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Abstract 



We derive a pair approximation (PA) for the NO+CO model with instanta- 
neous reactions. For both the triangular and square lattices, the PA, derived 
here using a simpler approach, yields a phase diagram with an active state 
for CO-fractions y in the interval y\ < y < y^, with a continuous (discontin- 
uous) phase transition to a poisoned state at y\ {yz). This is in qualitative 
agreement with simulation for the triangular lattice, where our theory gives a 
rather accurate prediction for yi . To obtain the correct phase diagram for the 
square lattice, i.e., no active stationary state, we reformulate the PA using 
sublattices. The (formerly) active regime is then replaced by a poisoned state 
with broken symmetry (unequal sublattice coverages), as observed recently 
by Kortliike et al. [Chem. Phys. Lett. 275, 85 (1997)]. In contrast with 
their approach, in which the active state persists, although reduced in extent, 
we report here the first qualitatively correct theory of the NO+CO model on 
the square lattice. Surface diffusion of nitrogen can lead to an active state in 
this case. In one dimension, the PA predicts that diffusion is required for the 
existence of an active state. 
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I. INTRODUCTION 



Following the introduction by Ziff, Gulari, and Barshad of a simple lattice model (ZGB) 
for the kinetics of the reaction CO + 1/2 2 ^ C0 2 on a catalytic surface, the study of sur- 
face reaction models has attracted increasing attention in nonequilibrium statistical physics 
IJ. Motivated by possible applications as well as intrinsic interest, the phase diagrams of 
a wide variety of models have been investigated in simulations and approximate, mean- 
field-like analyses. A typical feature is the existence of one or more absorbing states, i.e., 
configurations from which the system cannot escape 0. Continuous phase transitions to an 
absorbing state fall generically in the class of directed percolation While this aspect is 

highly universal, other details of the phase diagram depend on very specific model-dependent 
properties such as steric or geometric effects, the possibility of non-reactive desorption, dif- 
fusion of, and interactions among, adsorbed species. Applied to surface reaction models, 
mean-field theories, particularly at the two-site or pair level, often provide reasonable qual- 
itative predictions for the phase diagram. 

In the present work we derive pair mean-field approximations for one of the more com- 
plicated surface reaction models, that of NO+CO ||. We briefly review its main features, 
deferring a precise definition to Sec. II. The catalytic surface (i.e., one of the platinum- 
group metals), is modelled by a regular lattice (typically square or triangular) of equivalent 
adsorption sites. This surface is exposed to a reservoir of CO and NO at fixed concentra- 
tions. While CO needs but a single vacant site, NO requires a nearest-neighbor pair of sites 
to adsorb. (We note that a more realistic model, for example of the NO+CO reaction on 
Pt(100), would permit NO to adsorb at a single site; a vacant neighbor is required for disso- 
ciation JF].) The fundamental control parameter of the model is y, the probability that the 
next molecule arriving at the surface will be CO. One may also introduce nearest-neighbor 
hopping rates for the various adsorbed species. Nearest-neighbor CO-0 and N-N pairs are 
highly reactive: if any form (by adsorption or diffusion), they are eliminated before anything 
else happens, and the products (CO2 and N 2 , respectively) desorb immediately. 

On the triangular lattice (i.e., coordination number six), the phase diagram of the 
NO+CO model resembles that of the ZGB model: there is a reactive window for yi < y < y 2 , 
with a continuous phase transition at y± ~ 0.17, and a discontinuous transition at y 2 ~ 0.35 
PJ5|-|iT| . For y values outside the reactive window, the system eventually falls into an absorb- 
ing or "poisoned" configuration, devoid of vacant sites (the number of such configurations 
grows exponentially with system size). For y < y\ the final configuration consists predom- 
inantly of O, with an appreciable fraction of N. (A special case is y — 0, corresponding to 
a kind of random sequential adsorption (RSA) [|ll|] of dimers, with partial reaction. The 
final state consists of O and N atoms, with isolated vacancies interspersed.) For y > y 2 , CO 
takes over the role played by O in the small-?/ case. Diffusion of N, CO, and/or O shifts the 
transition points to some extent, but does not modify the phase diagram in any fundamental 
way. 

On the square lattice, the picture is radically different, there being (without diffusion of 
N), no active stationary state, whatever the value of y P,|l"2"|-|Ilf| . This observation, based 
on Monte Carlo simulations, was explained by Brosilow and Ziff (BZ), who argued that the 
active state is unstable to the filling of one sublattice with N atoms ||. Once this occurs, 
further adsorption of dimers is blocked, and the remaining vacancies are filled in with CO, 
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yielding a poisoned configuration. While BZ cast their argument in terms of global sublattice 
filling, in practice the system poisons through the growth of local "antiferromagnetic" do- 
mains, i.e., patches having one or another sublattice filled with N. Since the dynamics stops 
once the domains fill the system, relaxation to a globally-ordered state is not possible. For 
the same reason, no sharp transitions in coverages are observed as y is varied |^|.|12|.|13|.||.||.|T5[ . 
(It may be overstating the case to say that currently available results rule out any phase 
transition in the square lattice. Simulations show the coverages changing rapidly over a 
narrow range of y, but without discontinuities in the coverages or their slopes.) That the 
N-sublattice instability is responsible for destroying the active state (despite the absence of 
global sublattice order), is indicated by the observations that (1) an active state exists in 
the triangular lattice (which does not admit a decomposition into two sublattices), and (2) 
that diffusion of N (but not of O, or of CO) restores the possibility of an active state ||I5| , |I6| . 

Several theories of the NO+CO model have been proposed. Truncated at the 1-site level, 
the hierarchy of equations governing the cluster probabilities yields a reasonable estimate for 
y 2 on the triangular lattice ||||, but places the continuous transition at y\ = 0. Cortes et al. 
derived a pair approximation for the NO+CO model including CO-desorption, using finite 
reaction rates [F/J]. Kortliike, Kuzovkov and von Niessen (KKN) derived a very accurate 
prediction for yi on the triangular lattice using a two-site cluster approximation |18| ], 

None of the theories mentioned so far gives the phase diagram correctly for the square 
lattice: all predict an active state over some range of y. KKN made the fundamental 
observation that in this case, one must allow different concentrations on the two sublattices, 
to have any hope of capturing the instability identified by BZ |18| . They devised a two-site 
cluster approximation incorporating sublattices, and obtained an active state of reduced 
extent, and a third transition point inside the (mainly) CO-poisoned phase. In other words, 
the theory of KKN, while representing an improvement on theories ignoring sublattices, 
remains at variance with simulation (and the BZ argument), in allowing an active state, and 
predicts a third, unobserved transition. 

In this work we formulate a pair approximation (PA) for the NO+CO model on the 
square and triangular lattices, as well as in one dimension; we retain the instantaneous 
reactions generally used in simulations. (While a theory employing a finite reaction rate k 
is simpler algebraically, its k — > oo limit is not equivalent to a theory with instantaneous 
reactions |19fl .) The standard PA predicts an active state for < y\ < y < y 2 ; this is 
qualitatively correct for the triangular lattice, wrong for the square lattice. We obtain the 
correct phase diagram in the latter case from a pair approximation incorporating sublattices 
(PAS); the regime exhibiting activity in the PA now poisons via the sublattice instability. 
Our results for the coverages are in good accord with simulation, but due to the assumed 
homogeneity within each sublattice, phase transitions between different kinds of absorbing 
states persist at y\ and y 2 , and the order parameter (the difference in sublattice coverages), 
takes a nonzero value for yi < y < y 2 . We find that diffusion of N atoms lifts the instability, 
permitting an active state. But in our theory, the entire range yi < y < y 2 becomes active 
once the diffusion rate exceeds a critical value. In simulations the "active window" 
opens gradually as D n is increased |T5|JT^| . 

We have devised a simplified approach to deriving cluster mean-field equations. The 
method, which we illustrate with a simple example, proves particularly useful in the case of 
the NO+CO model, which allows eight kinds of nearest-neighbor pairs, and up to twenty- 
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two transitions among them. The remainder of this paper is organized as follows. Section 
II contains a definition of the model, including the reaction and diffusion steps. The PA 
method is described in Sec. Ill, with examples of its application to the contact process and 
the NO+CO reaction given in Appendices A and B, respectively. Our results are presented 
in Sec. IV, and a brief discussion follows in Sec. V. 



II. MODEL 

The NO+CO surface reaction model follows the Langmuir-Hinshelwood mechanism, in 



which both reacting species must be adsorbed on the substrate | 20fl . The steps below char- 
acterize the model ||: 

(a) CO (9) + V -+ CO a 

(b) NO (s) + 2V - O a + N a 

(c) 2N a -> N 2(9) + 2V 

(d) CO a + O a -> CO %) + 2V 

where A a indicates an adsorbed species, A( g ) a molecule in the gas phase, and V a vacant site. 
Steps (a) and (b) represent adsorption of carbon monoxide, and of nitric oxide, respectively. 
In step (c), two nearest-neighbor nitrogen atoms combine to form N 2 ( 9 ), and in step (d), an 
oxygen atom reacts with a carbon monoxide molecule to form C0 2 . N 2 ( ff ) and C0 2 ( 9 ) desorb 
from the surface immediately. In this work we assume complete dissociation of NO. (Various 
aspects of incomplete dissociation are considered in Refs. pl)|-|T7|j.) As noted above, reactions 
are assumed to occur instantaneously: nearest-neighbor CO-0 and N-N pairs cannot reside 
on the lattice. 

We now define the Markov process associated with the above set of reactions. On a 
lattice £ comprising N sites, the state-space of the process is the set of configurations 
{a} = {<Tj} ig £, where the site variable <7j takes values V, N, C or O in case site i is vacant, 
or occupied by N, CO, or O, respectively. One trial (or sample path) of the process consists 
of a sequence of configurations {cr} , {cr}^. A transition between configurations {cr} n 
and {<r} n+ i is generated via the following steps: 

(1) Choose the identity of the next arriving molecule: CO with probability y, NO with 
probability 1 — y. 

(2) In case of CO, choose a site x; in case of NO, choose a nearest-neighbor pair (x, y). If 
x (and/or y, in the case of NO) is occupied in {a} n , then {a} n+ i = {cr}„, i.e., the 
configuration does not change. Otherwise, let {a'} be the configuration obtained by 
placing CO at x, or in the NO case, N at x and O at y. 

(3) If {a'} contains no N-N or CO-0 nearest-neighbor pairs, then {er} n +i — {&'}■ If W} 
does contain such pairs, they will react. Specifically, if the newly-arrived CO has m 
neighbors in state O, then one of these (chosen at random if m > 1), as well as the CO, is 
removed from {a 1 } to give {cr} n+ i. In the case of NO deposition, the analogous procedure 
is applied to the newly-arrived N atom (if it has one or more neighbors N), and to the O 
atom (should it have any CO neighbors), to generate {a} n+ i- 

We associate with configuration {cx} n a "time" t = n/N. (This adds nothing to the pro- 
cess; it is convenient, nonetheless, to define a time unit comprising one attempted transition, 
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on average, per lattice site. In simulations, it is often more efficient to choose the first site 
x for the adsorption step from a list of currently vacant sites. But in this case, the time 
increment associated with the transition from {cr} n to {cr} n+ i is At = 1/V n where V n is the 
number of vacant sites in {cx} n .) 

To our knowledge, all of the simulation studies of the NO+CO cited herein treat the 
Markov process defined above. (The const ant- coverage studies of BZ clearly follow a different 
procedure, but the stationary properties of the two processes should converge in the large- 
size limit [§].) In the case of diffusion, however, various definitions have been employed. 
Here, we implement diffusion of (for example, N atoms) by modifying steps (1) — (3) as 
follows. Prior to (1), we impose 

(0) Choose the process: diffusion with probability D/(l + D), adsorption with probability 
1/(1 + D). In the latter case, proceed to steps (1) — (3) as above. In the former, perform 
instead: 

(1') Choose a site x at random. If x is not occupied by N, then {a} n+ i = {cr} n . Otherwise, 
choose a neighbor y of x at random. If y is occupied, {a} n+ i = {cr} n . Otherwise, let {a'} 
be {o~} n with x vacant and y occupied by N. 

(2') If {a'} is free of N-N pairs, {o"} n +i = {&'}■ Otherwise, choose a reacting pair, as in 
step (3) above, to generate {cr} n+ i. 

This procedure is not calculated to optimize computational efficiency, but rather to provide a 
meaning for the parameter D that is independent of the configuration. The rate of hopping 
attempts of an adsorbed N atom is D/(l + D). Diffusion processes for other species are 
defined analogously. 



III. PAIR APPROXIMATION 

Before considering the NO+CO model in detail, we explain a simplified method for 
deriving the PA equations. These equations govern the evolution of the probabilities P(ij, t), 
that a randomly chosen nearest-neighbor pair of sites, x and x', say, are in states i and 



j. In most previous treatments, (see, for example Ref. pi]), the rates of change of the 
P(ij) are derived by enumerating the changes in the number of nearest-neighbor pairs in 
a neighborhood of sites centered on, and including, a central pair. For example, in a one- 
dimensional system with nearest-neighbor interactions and two states, '0' and '1', per site, a 
transition of the form (0011) — > (0111) occurs at rate P(0011)w(001 — > 011). Counting the 
changes in the central pair and the periphery, we see that in this process, one (11) pair is 
created and one (00) pair destroyed. The simplification comes from ignoring changes outside 
the central pair, and regarding the above process as one in which (01) — > (11). Since P(ij) 
is the probability for any nearest-neighbor pair, following the changes at a particular pair 
(e.g., the central one), is sufficient. This results in a significant reduction in bookkeeping, 
particularly in two or more dimensions, and for processes (such as the NO+CO model) in 
which a fairly large number of peripheral sites can influence the transition probabilities. We 
illustrate the method by applying it to the contact process in Appendix A. 

In the table below we list the allowed (-y/) and forbidden (0) nearest-neighbor pairs in 
the NO+CO model. (Entries below the main diagonal are redundant.) 
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Next we require the set of transitions between pairs. In the table below, we assign 
arbitrary labels to the allowed transitions, and leave the remaining fields blank. For the 
square lattice only processes 1-20 are pertinent, 21 and 22 being possible only on the 
triangular lattice. In one dimension, transition 4 is also excluded. Diffusion alters the rates, 
but not the set of possible transitions, the only exception being when we consider sublattices. 
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The heart of the calculation lies in deriving expressions for the transition rates R±, i?22- 
Once these are in hand, we can write the equations for the pair probabilities by noting that 
each transition acts as a source for one pair (i.e., in the first column of the table), and a 
loss term for another (listed in the top row). Denoting the pair probabilities by (ij), where 
i and j can be V, N, C, or O, we have, for example, that 



d(VN) 
dt 



-R5 + Re + R7 + R21 — R\ — R\a — R17 ~~ R 



22- 



and 



d(VV) 
dt 



— 2[Ri + i?2 + -R3 + R4 ~~ R5 — Rs ~ R\\ 



R 



16j, 



(1) 



(2) 



the overall factor of two arising, in the last expression, due to contributions in which the 
molecules (always of different species, in this case) occur in the opposite order. Here it is 
important to emphasize that for i 7^ j, (ij) represents the probability to find a site x in state 
i, and its neighbor y in state j; (ji), which represents the reversed situation, is of course 
equal to (ij) by symmetry. Thus the normalization condition reads 



(VV) + (CC) + (00) + 2[(VN) + (VC) + (VO) + (NC) + (NO)] = 1. 
We introduce a similar notation for site probabilities: 

(V) = (VV) + (VN) + (VC) + (VO), 



(3) 



(4) 



(N) = (NV) + (NC) + (NO), 



(5) 
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(O = (CV) + (CN) + (CC), (6) 



and 

(O) = (OV) + (ON) + (00). (7) 

Another useful piece of notation represents the probability of having site x in state i and its 
neighbor not in state j by (if). For example, 

(VM) = (VV) + (VO) + (VC), (8) 

and 

(C(p) = (CV) + (CN). (9) 

Finally, when employing sublattices A and B, we use (i)a to denote the site probability 
of i in sublattice A (similarly for (i)s), and (i])a to denote the probability of finding a site 
in state i in the A sublattice, and its neighbor (in B), in state j. (if) a is defined analogously. 
In a sublattice calculation we have thirteen different pair probabilities, since for i ^ j we 
must distinguish (ij)A and (ij)s- (By definition, (ij)A = (j^)b-) Diffusion introduces one 
further transition beyond those enumerated above: if species i can perform nearest-neighbor 
hopping, the transition (iV)A — ► (iV)s becomes possible. 

Each (nondiffusive) event involves the arrival of a molecule, either CO or NO, at the 
surface. Since the next arriving molecule is CO with probability y, rates for processes 
involving the arrival of CO carry a factor of y. In processes involving the arrival of NO, we 
require the probability that the next event involves N arriving at a certain site x, and, of 
course, O arriving at a neighbor, y. In a lattice with coordination number z, this probability 
is y = (1 - y)/z. 

The expressions for the various rates are, in general, quite complicated, and we shall not 



list all of them here. Examples of their derivation are given in Appendix B The PA 
equations are integrated numerically, using a fourth-order Runge-Kutta scheme [23], starting 
from an empty lattice. 



IV. RESULTS 
A. One dimension 

Applied to the NO+CO model on a line, the PA predicts no active steady state in the 
absence of diffusion (see Fig. 1). The dependence of the coverages on y is qualitatively 
similar to that found in simulations of the square lattice P,|T5|. The vacancy fraction is 
nonzero only for y = 0, where we find (V) = 0.1623, (O) = 0.5013, and (N) = 0.3364 in 
the stationary state. In fact, the final vacancy concentration should be the same as in one- 



dimensional dimer RSA (without reaction), i.e., (V) = e = 0.13534... P^| . This is because 
a N-N reaction always yields the configuration OVVO. The vacancy pair is subsequently 
filled in, so sum of the final coverages, (N) + (O), is identical to the final O coverage 



in RSA of O2 p5 |. While the PA is exact for dimer RSA in one-dimension (due to a 
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shielding property ]TT[] ) , in the present case the approximation cannot deal adequately with 
the reactions. The configuration VNV, for example, is impossible in the pure-NO process 
(each N has at least one O neighbor), but is assigned a nonzero probability in the PA. We 
obtain better results for the final coverages from a three-site approximation: (V) = 0.1486, 
(O) = 0.5001, and (N) = 0.3505. (The KKN method yields a further slight improvement: 
(V) = 0.1453, (0) = 0.5095, (N) = 0.3452; simulations yield (V) = 0.1353, (O) = 0.5066, 
and (N) = 0.3581 J26J.) 

Next we consider the effect of a nonzero diffusion rate, D^, of N a atoms. For D^ > 
D C N = 4.38, we find a reactive window for y 1 < y < y 2 , with a continuous (discontinuous) 
transition at yi (y 2 ). Fig. 2 shows the coverages for D N = 10.0. In Fig. 3 we plot A = y 2 — y\ 
as a function of D^. With increasing Dn, yi tends to zero and y 2 to 0.2. Close to D c Nl A 
is described, approximately, by A ~ (Dn — D N ) 067 . F° r diffusion of CO (but none of the 
other species), we find an active state for Deo > D co = 7.19. The width, which is very 
small, grows linearly with Deo m the neighborhood of the critical value. Both y\ and y 2 
shift to higher values with increasing Deo- Finally, for diffusion of O atoms (exclusively), 
we find Dq = 0.75. The window width follows A « (Do — -D^) ' 24 in the vicinity of D C Q . 
The continuous transition (y = y\) always occurs near y — 0; for large values of Do, yi 
approaches a limiting value of 0.27. We note that for y — 0, even very small values of Dq 
cause a drastic reduction in the fraction of vacant sites: when Dq = 0.05, for example, 
(V)< 10~ 4 ; for Dq = 16.0, (V) < 10~ 6 . Similar behavior is observed in simulations on the 
triangular lattice [p? . 



B. Triangular lattice 



The coverages predicted by our method for the NO+CO model (without diffusion) on the 
triangular lattice are compared against simulation results || in Fig. 4. There is an active 
steady state between the continuous transition at y\ = 0.040(1), and y 2 = 0.363(1), which 
marks a discontinuous transition. (We note that the latter is the result for an initially empty 
lattice. Commencing with a finite CO-coverage will in general yield a smaller value for y 2 .) 
The table below compares our PA results for the transition points with those predicted by 
the site approximation (SA) ||||, and the KKN method |18|], and found in simulations [§]. 





SA 


KKN 


PA 


Simulation 


yi 


0.0 


0.152(1) 


0.040 


0.1725(25) 


V2 


0.3877 


0.393(1) 


0.363 


0.3514(1) 



While our result for y 2 is in good agreement with simulation, we obtain a poor estimate for 
yi. The latter is a consequence of neglecting explicit correlations beyond nearest neighbors 
in the pair approximation, and, perhaps, of ignoring certain nearest-neighbor pair factors 
in reckoning cluster probabilities (see Appendix B). For y = we find a final poisoned 
state characterized by (O) = 0.8538, (N) = 0.1462, and (V) ~ 2 x 10~ 6 . Simulations yield 
(O) = 0.7443, (AO = 0.1656, and (V) = 0.0901 0. Clearly the PA does not give an 
accurate description of this RSA process with partial reaction, for either the triangular or 
square lattices (see below). (This is in contrast to pure dimer deposition, for which the PA 
does reasonably well ||21| .) Despite these discrepancies, the PA coverages are generally in 
good agreement with simulation. 
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C. Square lattice 



The PA prediction for the phase diagram is qualitatively similar to that found for the 
triangular lattice. The continuous transition from a predominantly O-poisoned state to an 
active state occurs at y\ = 0.111, and the discontinuous transition falls at y 2 = 0.2981. (The 
latter, again, is determined using an initially empty lattice.) Cortes et al. obtained y\ ~ 0.09 
and y 2 ~ 0.35 in this case, showing that the large-reaction-rate limit of a calculation using 
finite rates yields results comparable, but not identical to, one employing instantaneous 
reactions. (While the comparison is academic in the present instance, there being no active 
state on the square lattice, it is of interest to gauge the agreement between the two methods.) 
For y = we obtain (O) = 0.6551, (N) = 0.0198, while simulations yield (O) = 0.6495 and 
(iV) = 0.2416 |27|. 

Since the PA predicts no active state in one dimension (in the absence of diffusion), it 
is interesting to check whether removing reaction 4, which is impossible in one dimension, 
changes the phase diagram. We find that deleting this process has a minimal effect on the 
PA prediction for the square lattice. The cause for the dramatic difference between the 
linear and square lattice results must be sought elsewhere. (It is worth remarking that the 
PA similarly yields only absorbing states for the one-dimensional ZGB model.) 

Introducing sublattices in the PA calculation (yielding what we call the PAS theory), 
changes the result drastically. We find that the active state predicted by the PA is unstable 
to the formation of N-rich and N-poor sublattices; in the process, the vacancy density falls to 
zero, and the active state vanishes. As detailed in Appendix B, we employ an extended set 
of variables, i.e., for % 7^ j, = V, N, C, or O), probabilities {ij)a an d (ij)B, representing 
species i in the A or the B sublattice. (Naturally, the site coverages (N)a and (N) B , etc., may 
also differ.) We begin the calculation, as before, with an empty lattice; the pair equations 
reach the same steady solutions as in the simple PA. For yi < y < y 2 , we then probe the 
stability of the active state by transferring a small amount, A (AT) = (N)a — (N)b = 10~ 3 , 
of N from one sublattice to the other, and study the response to this small perturbation. 
We find, in all cases, that A(N) grows, and that (V) decreases, finally becoming zero, that 
is, the PAS equations reach an absorbing state. Essentially identical results are obtained if 
we start from a slightly asymmetric initial condition, i.e., with a small N-coverage on one of 
the sublattices, as was done in Ref. |TJ|. (We note in passing that introducing sublattices 
in the site approximation has no effect on the results.) 

The PAS, then, represents the first theoretical approach giving a phase diagram in accord 
with simulation, for the NO+CO model on the square lattice. But it retains some of the 
undesirable features of the PA. Fig. 5 shows that for y < 0.1, and again for y > 0.5, the 
PAS coverages are in good agreement with simulation. In these regimes, of course, the PA 
and PAS are identical (they only differ on the interval where the PA predicts an active 
state). For 0.1 < y < 0.5 there are substantial differences between theory and simulation, 
associated with the continued appearance of phase transitions at y\ and yi in the PAS. For 
y\ < y < y 2 the PAS equations exhibit spontaneous symmetry breaking [[Tj|] associated with 
global "antiferromagnetic" order, A(N) 7^ 0, (see Fig. 6), which, as we have remarked, is 
not seen in simulations. 

A final point concerns the effect of N diffusion on the phase diagram. Since this process 
tends to equalize the sublattice coverages, it is reasonable to expect an active state for 
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sufficiently large values of D^. This is indeed observed in simulations [fT5p8f , where the 
range A of y values supporting an active state grows steadily with D N . (The PA, as we 
have noted, yields a similar result in one dimension.) When we include N diffusion in the 
PAS calculation, we observe no active state for D N < D C N ~ 0.023, but for larger diffusion 
rates the entire interval between y\ and 2/2 becomes active at once. 

V. DISCUSSION 

We have formulated a pair approximation (PA) for the NO+CO surface reaction model 
with instantaneous reactions, using a simplified derivation. The PA gives quite reason- 
able predictions for coverages on the triangular lattice, but a surprisingly low value for the 
continuous transition point, y±. The pair approximation with sublattices (PAS) gives a qual- 
itatively correct phase diagram for the square lattice, but certain anomalous features of the 
simple PA persist, notably, singular behavior of the coverages at the transition points yi and 

Our study shows that the PA can furnish reliable qualitative, and in some instances 
quantitative predictions for reaction models, provided it includes a mechanism for realizing 
any symmetry-breaking tendency inherent in the model. The same condition applies, of 
course, to mean-field or cluster approximations for equilibrium models. In general, it is 
asking too much of such theories to provide quantitatively reliable phase boundaries; the 
PA may nonetheless yield some insight into the overall shape of the phase diagram. More 
accurate theories are typically based on the hierarchy of n-point functions. 

The theory of KKN yields a remarkably accurate prediction for yi on the triangular 
lattice. This indicates that the effect of long-range correlations is reasonably well-represented 
in their theory. Surprisingly, the PA yields a better prediction for 7/2. It remains to develop 
a method that combines the advantages of the two approaches. 
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Appendix A: Simplified derivation of pair equations 



Here we illustrate our method by applying it to the contact process (CP) ||29|| . The CP 
is a Markov process defined on a <i-dimensional cubic lattice. Each site is either vacant (0) 
or occupied (1). The transition rates at any site are w(l —>■ 0) = 1 (independent of the 
neighbors) and w(0 — > 1) = Xn/2d, where < n < 2d is the number of neighboring sites 
in state 1. (Note that all sites is an absorbing configuration.) Thus in one dimension 
we have w(101 -> 111) = A, w(100 -> 110) = w(001 -> Oil) = A/2, and w(000 -> 
010) = 0. We enumerate below the transitions, associated rates (e.g., the transition rate 
times the probability of the initial configuration), and overall changes in the number of l's 
and 11 nearest-neighbor pairs, in the one- dimensional CP. We use (1) to denote the density 
of l's, (10) for the probability of a nearest-neighbor — 1 pair, etc. (By normalization, 
(00) + 2(01) + (11) = 1. Note that the second and fourth entries carry a factor of 2 to 
account for mirror-image events.) 



Process 


Rate 


AiVi 


AiVn 


101 - 


-> 111 


A(01) 2 /(0) 


+1 


+2 


100 - 


-> 110 


2A(01)(00)/(0) 


+1 


+1 


111 - 


-»■ 101 


(n)7(i) 


-1 


-2 


110 - 


-> 100 


2(11)(01)/(1) 


-1 


-1 


010 - 


-> 000 


(oi)7(i) 


-1 






Collecting results, one finds 



d(l) 
dt 



A[(l)-(11)]- (1), 



and 



rf(ll) 
dt 



A 



'D-(ll^ 



(0) 



[1 - (11)] " 2(11). 



(10) 



(11) 



In this approximation, an active stationary solution (one with (1) > 0), exists only for 
A > A c = 2. The above calculation is simple enough in one dimension, but becomes more 
complicated for higher d. We illustrate our simpler alternative method below on the d- 
dimensional CP. Only transitions at the central pair are enumerated. The rates involving 
creation receive independent contributions from within the pair (if it is of type 01) and from 
the 2d — 1 neighbors outside. (This independence is of course an approximation intrinsic to 
the PA.) 



Transition 


Rate 


11 -> 01 


Ri = (11) 


01 -> 11 


R 2 = A(01) [1 + (2d - l)(01)/(0)] /2d 


10 -> 00 


R 3 = (10) 


00 -> 10 


R A = (2d- l)A(00)(01)/2d(0) 



We find d(ll)/dt = 2[R 2 — R±], which reduces to Eq. (|TTJ) for d — 1, when we note that 
(01) = (1) - (11). Noting that (1) = (11) + (10), we immediately recover Eq. (|10|) for any 
d. (Analysis of the stationary solutions shows that A c = 2d/(2d — 1) in d dimensions.) 
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Appendix B: Rates for the NO+CO model 



In this Appendix we present several examples of the pair approximation (PA) transition 
rates for the NO+CO model, first on the square, and then on the triangular lattice. We begin 
with transition 15, VC — > NC, on the square lattice. Fig. Bl shows one of three equivalent 
configurations needed to realize this process. The rate (per bond of the lattice) carries the 
factors (VC) (probability of the initial state) and y (probability of NO arriving with N 
falling at the central vacant site). There are three possible locations (neighbors of V) at 
which O might adsorb; each is vacant, in the PA, with probability (VV)/(V). (Whether the 
O atom reacts or not is unimportant in this instance.) If either of the remaining neighbors 
of V harbors an N atom, the newly- arrived N will react. (Recall that we are considering 
infinite reaction rates in this work.) Hence these two sites must be free of N for the desired 
transition to occur, implying a factor of [(Vlfi) / (V)] 2 . Thus 

JVC)(VV)(V^) 2 

Rl5 = 3y (yf ' (12) 

Now suppose that N atoms can hop to neighboring sites at rate D^- If one of the neighbors 
of V harbors an N atom, it can hop to the vacant site, and will remain there if the other 
neighbors are free of N. We must then add to the above expression the diffusive contribution 

3D N (VC)(VN)(V^) 2 

Dl5 - ~1 (vf ' (13) 

the factor of 1/4 reflecting the four possible directions of the attempted hopping move. (Note 
that we absorb a factor of 1 + into a rescaled time variable.) 

If we consider sublattices, each process splits in two. For process 15, the rate (including 
diffusion) for the case in which V lies in the A sublattice is readily seen to be 

p o~ (yCUVV)(VP) 2 A , SD N (VC) A (VN) A (V^) 2 A 

R^A = 3y ^ + ^ . (14) 

(For any i, Ri^s is found by interchanging A's and B's in Ri tA ) 

A somewhat more complicated transition is number 4, NC — > VV. It is contingent upon 
an NO landing parallel to the central NC pair, with the two N's adjacent. One of two 
equivalent initial configurations is shown in Fig. B2. The solid lines indicate the pairs 
included in reckoning the probability of the configuration. Notice that when a pair of 
neighboring sites, x and y, are both neighbors of sites in the central pair, we use the pair 
factor associated with the central pair [(VN)(VC) in this example], in preference to the factor 
between x and y [either (VN)(VV) or (VC)(VV)]. We apply this rule in all our calculations, 
to eliminate possible ambiguities. Consider the neighbors a and b of the vacant site above 
N. If neither of these bear N, the newly-arrived N will surely react with the N in the central 
pair. If either a or b (but not both) carry an N, the probability of the desired reaction is 
1/2; it is 1/3 in case both a and b bear N. The probabilities of these events — no N, one N, 
or two — given the vacant neighbor, are (V)l) 2 /(V) 2 , 2(Vft)(VN)/(V) 2 , and (VN) 2 /(V) 2 , 
respectively. An analogous consideration applies to the probability of reaction of the central- 
pair CO, depending upon the presence of CO at c and/or d. Combining these observations, 
we arrive at 
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n _ (NC)(VN)(VC) 
4 " V (N)(C)(V)* 



(VN)(V) + -(VNf 



(W)(v) + -(vcf 



(15) 



where we used (VN) + (VN) = (V) and (V<?) + (VC) = (V). For the sublattice calculation, 



R. 



4,A 



2y 



(NC) A (VN) B (VC) J 
(N) A (C) B (Vf A (V)\ 



{VN) B (V) B + -(VNf B 



(VG) A (V) A + -(VC)\ . (16) 



There is no diffusive contribution to this process. 

While several rates have more complicated expressions, all of the calculations in the 
square lattice follow the lines of those illustrated above. A new question of principle does 
arise in the triangular lattice, where it is not possible to include all the pair factors between 
the central pair of sites and their nearest-neighbors in the simple PA. In the triangular 
lattice, if x and y are nearest neighbors, they have two neighbors (sites d and h in Fig. B3), 
in common. In the PA, the probability of finding x, y and d in states i, j, and k, resp., may 
be written as (ij)(ik)/(i), (ij)(jk)/(j), or (ik)(kj)/(k), but not as (ij)(ik)(jk)/(i)(j)(k). 
(The latter family of expressions is not, in general, even normalized!) Our choice of which 
pair factors to include is shown in Fig. B3. Note that as we sum over all the possible states 
of the peripheral sites, symmetry under the interchange of x and y is restored. 

As an example consider process 19, VC — > CC (site x vacant and y occupied by CO in 
Fig. B3). A CO molecule must land at x and remain there, which implies that sites a, b and 
c must be free of O. (Sites d and h have no possibility of bearing O, as they are neighbors 
of a site occupied by CO. The states of sites e, f, and g are unimportant in this process.) 
Multiplying the independent factors associated with the events enumerated above, we find 
that the rate for this process is given by: 



R 



19 



y{vc) 



(VQf 



(17) 



A more complicated process is number 20, VO — > 00 (x vacant, y occupied by O in 
Fig. B3). An NO must fall with O at x, and N at one of the sites in {a,b,c,d,h}. In no 
case may any of the sites in {a,b,c} hold CO. (Being neighbors of an O, d and h are surely 
free of CO.) The rate is given by the expression: 



R 



20 



y{ ] \ ivyyivy) + W)\WY) + ioT\ivy) 



(18) 



The first term in brackets represents N falling at a, b, or c. The next is for N falling at 
d. The probability of d being vacant, given one vacant neighbor and one occupied by O, 
is (VV)/[{VV) + (VN) + (VO)] = (VV)/(V(P). The final term represents N falling at h. 
These examples illustrate the principles used in the calculation, the resulting expressions, 
needless to say, becoming quite involved in certain cases. 
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Figure Captions 

FIG. 1. Coverages versus y in one dimension in the absence of diffusion. Solid line: (C); 
dotted line: (O); dashed line: (N). 

FIG. 2. Coverages versus y in one dimension, for Dn = 10.0. Symbols as in Fig.l. 

FIG. 3. Width of the reactive window in one dimension, as a function of the diffusion rate 
Djy. The critical value D C N = 4.38. 

FIG. 4. Coverages in the triangular lattice. Solid lines: simulation (Ref. M), dashed lines: 
PA. 

FIG. 5. a) N coverage in the square lattice. Solid line: simulation (Ref. ||), dashed line: 
PAS. b) CO coverage, symbols as in a); c) O coverages, symbols as in a). 

FIG. 6. PAS sublattice order parameter, A(N). 
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